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1. Introduction 

It is a long known fact, that interaction of a system 
with electromagnetic radiation is affected not only by 
the system itself, but also by the environment in which 
it is situated. Placement of an atom or a molecule in 
a cavity leads to many effects not present in unbounded 
vacuum, which are the subject of studies within the 
area of so-called cavity QED - for a review see e.g. [1]. 
In particular, one of these phenomena is the modifica- 
tion of spontaneous emission rate, predicted by Purcell 
[2] and later verified in many experiments, e.g. [3-6]. 
This fact has considerable consequences for contempo- 
rary science and technology, therefore it is an impor- 
tant object of research. Because spontaneous emission 
is a quantum effect, for its modeling it is necessary 
to employ a formalism based on quantization of elec- 
tromagnetic field. A theoretical description suitable 
for the treatment of spontaneous emission in dielectric 
structures, followed in many papers, has been discussed 
by Glauber and Lewenstein [7]. 

Another fact, that has been recognized over a hundred 



years ago and described by Lord Rayleigh [8], is that 
periodically arranged medium has peculiar properties, 
with the most interesting one being the existence of 
band gaps - frequency ranges in which no propagating 
waves exist. In photonics, periodic dielectric materials 
with band gaps for electromagnetic waves (photonic 
band gaps) are called photonic crystals. The simplest 
of this kind of structures are one-dimensional photonic 
crystals, built of alternating layers with two different 
refractive indices and widths (and, in fact, can be con- 
sidered a particular type of planar multilayer waveg- 
uide). They are a particularly good subject of study, 
because one-dimensional multilayers can be described 
analytically and are significant from the practical point 
of view. As it has been shown, by proper choice of pa- 
rameters it is possible to obtain a band gap for all di- 
rections of propagation [9-11], but even without a full 
band gap, one-dimensional periodical lattice allows for 
a significant modification of emission of light and there- 
fore is useful for construction of devices such as light 
emitting diodes or lasers [6, 12, 13]. To understand 
their operation, and to be able to design them with 
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care for various details, it is important to properly de- 
scribe the impact of the structure on the process of 
emission of light. 

A perfect photonic crystal, which is strictly periodic 
structure, consists of an infinite number of elementary 
cells (i.e. groups of layers, which appear in the same 
sequence in every period) and can be modeled with the 
use of Floquet-Bloch theorem, fixing the form of modes 
of the structure. A recent example of this treatment is 
[14], in which spontaneous emission has been described 
with a classically derived formula. However, a practi- 
cal realization of a photonic crystal is inevitably of fi- 
nite size, thus, in a more realistic models, the structure 
must be treated rather as a planar multilayer waveg- 
uide. This approach makes it possible to account for 
defects (of width or refractive index) of the structure 
as well. Just as in a perfect photonic crystal, it is 
possible to use a classical description of spontaneous 
emission in a planar waveguide, e.g. [15]. One has to 
note, though, that what is actually calculated this way, 
is radiation of a classical electric dipole, and only com- 
parison with a result from quantum theory allows to 
establish its relation to spontaneous emission. Hence, 
such treatment is rather just a limited reconstruction 
of a specific formula, and not a relevant, general theory 
in any way. A proper, versatile description of sponta- 
neous emission has to be established in the quantum 
framework, like [7]. For a multilayer waveguide, a par- 
ticular case of which is a finite one-dimensional pho- 
tonic crystal, one cannot use the Floquet-Bloch theo- 
rem and has to choose modes of the structure in a dif- 
ferent way. A basis suitable for a single interface be- 
tween two media has been proposed by Carniglia and 
Mandel [16]. In their proposition, the set of modes is 
constructed starting with plane waves propagating to- 
wards the structure (sweeping through all frequencies 
and angles of incidence), thus it can be referred to as 
a set of so-called incoming modes. Each incident plane 
wave has to be accompanied by respective reflected and 
transmitted waves, and such triplet constitutes a field 
distribution in the structure. Modes of this form are 
particularly simple, form orthogonal [17] and complete 
[18] set. There also exists a modification of this model, 
in which instead of waves propagating towards the in- 
terface, waves propagating away from the interface are 
used [19], so-called outgoing modes. With this choice 
there is always only one emitted plane wave carrying 



power out in a chosen direction, and this set is con- 
sidered favorable for this kind of calculations, e.g. in 
[20]. For a multilayer structure, it is necessary to sup- 
plement these models with guided modes, calculated 
for the layer with the highest refractive index [21], and 
radiation modes can be redefined as combinations of 
the incoming or outgoing modes, for which no power 
flow in transverse direction occurs, what assures their 
orthogonality [22]. Field distributions of modes can be 
conveniently found with the help of translation matrix 
method [23, 24]. 

For all these kinds of modes, it is a common feature, 
that a kind of source of radiation exciting a mode is 
situated outside the structure, what seems not to be 
particularly suitable for description of emission from 
a system located inside the structure, because it leads 
to treatment of the structure as a whole block, instead 
of distinguishing particular locations inside, which may 
have significantly different properties. However, it is 
possible to find field distributions starting from a plane 
wave emerging in one of the layers. This approach, de- 
scribed or employed e.g. in [25, 26], requires to sum 
up all the reflected waves appearing in the layer, which 
resembles a text-book exercise of calculating transmis- 
sion through a FP resonator. The difference is that in 
the case of source inside the layer, reflections occur at 
both interfaces, thus there are twice as much waves to 
sum. Modes obtained with this method are indexed 
by wave vectors and polarizations, and naturally split 
into radiation, substrate and guided modes. This con- 
struction of modes is also the starting point of the ef- 
fective resonator model [27-31], which allows a very 
clear physical interpretation of obtained results. The 
model is not just limited to calculating distributions of 
electromagnetic field. Its main advantage is the defini- 
tion of a quantity called mode spectrum, which can be 
used to quickly and conveniently analyze modes and 
properties of the structure. Assuming, that the inter- 
faces between layers are perpendicular to z axis, mode 
spectrum of a specific layer is defined as 

WM- 1 l-\r L r R \ 2 

8tt 3 \l-r L r R exp(2ik z L z )\ 2 ' V ' 

where k is wave vector and e - polarization of the mode, 
tl and tr are reflection coefficients of the stacks of 
layers on both sides of the layer with the source (i.e. 
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subscript L stands for "left" and R - for "right", these 
coefficients are functions of polarization and wave vec- 
tor, but for clarity of the expression it is not explic- 
itly denoted), k z is the z component of the wave vec- 
tor and L z is the layer's width. Similar expressions 
have already appeared in the literature, but they have 
never been given such an important physical mean- 
ing. Mode spectrum is a quantity directly related to 
density of states and spontaneous emission rate, resem- 
bling quality factor of a cavity, which also takes into 
account interference between plane waves composing 
a mode. If the quality factor for a mode is high, then 
the value of mode spectrum is high when the inter- 
ference is constructive, or low, if the interference is 
destructive. In particular, if mode spectrum is equal 
to zero, this means, that the field of the mode is com- 
pletely extinguished in the cavity. For modes, for which 
the quality factor is low, the value of mode spectrum 
is close to free space value 

In this paper, we derive the expression for the spon- 
taneous emission rate in a one-dimensional photonic 
crystal with arbitrary defects of width or refractive 
index of any layers, split into explicit contributions 
from different kinds of modes, i.e. radiation, substrate 
and guided modes, obtained with the help of the ef- 
fective resonator model. It is shown, that mode spec- 
trum characterizes spontaneous emission rates inher- 
ent for each layer. It allows for an easy investigation 
of properties of the structure, in particular containing 
defects, either introduced intentionally [32] or random 
fabrication impurities, which have negative influence 
on the structure [33]. Calculations conducted with ef- 
fective resonator model lead to description of spon- 
taneous emission which is particularly convenient for 
research on impact of the structure on the process or 
dealing with practical tasks, where a proper design of 
the structure is necessary. Obtained results can be 
used for a detailed study on how a multilayer environ- 
ment affects spontaneous emission decay, first of all to 
identify the modes which contribute the most in par- 
ticular layers, what is not straightforward with other 
construction of modes. Effective resonator model is 
completely analytic, thus, in principle, it allows to ob- 
tain analytic expressions describing various aspects of 
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Figure 1. Sketched structure of dielectric multilayer waveg- 
uide. 



the phenomenon or could be used as a basis for semi- 
analytic or approximate calculations in more compli- 
cated structures. The paper is organized as follows. In 
Sec. 2, we define the considered structure, introduce di- 
mensionless parameters, particularly suitable for anal- 
ysis of its properties, and briefly review a few proper- 
ties of modes in the effective resonator model. Next, 
in Sec. 3, we describe the quantum formalism adopted 
for calculations, which we employ in Sec. 4 to obtain 
the expression for spontaneous emission rate. Contri- 
butions from different kinds of modes are the subject 
of calculations in Sec. 5. We provide a few exemplary 
results in Sec. 6 and Sec. 7 summarizes the paper. 

2. Considered structure and 
adopted notation 

We consider a structure of finite one-dimensional pho- 
tonic crystal, which is a particular case of a dielectric 
multilayer waveguide, sketched in Fig. 1. We define 
the z axis perpendicular to interfaces between dielec- 
tric layers, which we assume to be made of lossless, 
isotropic dielectrics. The refractive index of the struc- 
ture is then given by 

n ( r ) = X)xi(*)rau), ( 3 ) 

where j indexes layers, nu-\ is the refractive index of 
the jth layer, and Xj( z ) the characteristic function of 
the jth. layer: Xj( z ) = 1 for « in the jth layer and 
Xj( z ) = for z from outside of it. Index j runs from 
-Nl < to Nr > and we put z = at the left 
boundary of the layer for which we choose j = 0. We 
treat the media outside as layers with j = — Nl and 
j = Nr, but additionally assume, that ri(R) < ri(L) 
(in subscripts and superscripts we use abbreviations 
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R = Nr and L = —Nl) - this is just the matter of z 
axis orientation, therefore it does not affect generality 
of the model or presented calculations. We denote the 
width of the jth layer by Lg , using for the j ' — layer 
L z = Lf\ and assume that layers extend to infinity 
in x and y directions. For one-dimensional photonic 
crystal of our concern 



Hi) 



ni for even j, 
n2 for odd j, 



and 



. (j) | Li for even j, 



L-2 for odd j, 



(4a) 



(4b) 



except for defected layers and the surrounding media, 
for which the general notation n(j) and has to be 
kept. 

Modes of the structure are indexed with wave vector 
k and polarization e = TE,TM. To each of them 
there are bound electric and magnetic field distribu- 
tions, fke{r) and f" e (r). Different layers are equipped 
with different sets of modes, in the paper we always 
refer to modes of the layer for which j = 0, thus we 
obtain results relevant for this particular layer (but 
chosen arbitrarily). Construction of these modes is dis- 
cussed with details in [28-31], here we wish to briefly 
review only a few of their properties. 
Electric and magnetic fields radiated by a source in- 
side the j = layer can be written in terms of field 
distributions f^ e and f^ e : 

E{r,t) = t £ke(t)f k e(r) (5a) 
= V / d 3 fc£ fce (t)/ fce (r) 

e J RSM 

+E/ d2fc n E &„«(*)/*.«(»•). 

e aSGM(fe||,e) 

H(r,t)=£ B fce (*)/£(r), (5b) 

where £ke(t) and Bke(t) are time-dependent ampli- 
tudes, RSM denotes subset of radiation and substrate 
modes, and GM(fe||,e) - discrete subset of guided 
modes with given fey = \Jk% + k% and polarization e. 
For guided modes the integration is performed over 
components k x and k y , while the whole wave vector of 



the ath guided mode is k a = k x e x + k y e y + k Zta e z . 
Decomposition of modes into radiation, substrate and 
guided modes is based on a standard criterion of total 
reflection, we precisely define each of these subsets in 
Sec. 5. To improve clarity, we use the "sum-integral" 
symbol, which combines integration over components 
of wave vector and summation over discrete guided 
modes - expression (5a) can be considered its defini- 
tion. 

Every mode can be identified with its electric field dis- 
tribution, because magnetic field distributions can be 
obtained from the relation: 



where 



fH 

J fee 



V X fke 



c 

"(0) 



'k 2 



(6) 



(7) 



is the angular frequency of the mode. The last expres- 
sion is the well known dispersion relation for a plane 
wave, which is a condition the plane wave must sat- 
isfy to be a valid mathematical solution of Maxwell 
equations, in which LJk and components of k can be in 
general complex numbers. Further restrictions can be 
settled by requirement, that the field must be physical, 
i.e. must not "explode" in infinity, meaning that the ex- 
ponents exp(ife • r) exp(— icjfet) should either oscillate 
or fade. An obvious notion is that u)k £ K, because 
of the time uniformity, similarly k x ,k y G R, because 
the considered structure is uniform in x and y direc- 
tions, but spatial dependence in z direction is more 
complicated. The field is always physical if k z G R, 
but in general it can be as well for some imaginary 
values k z G iR, if the refractive index of the layer is 
sufficiently small (the field has to be oscillatory in an- 
other layer, to connect two evanescent waves fading in 
opposite directions). However, in this paper we wish 
to concentrate only on modes with real wave vectors 
k G R 3 , i.e. our considerations are limited to layers 
with sufficiently high refractive index. The reason is 
that for a single non-uniform plane wave with imagi- 
nary k z there is no flow of energy in z direction, thus, 
a natural choice of waves on an interface between two 
layers does not contain a reflected wave. As our defi- 
nition of mode spectrum is based on reflection coeffi- 
cients, it cannot be naturally applied for non-uniform 
plane waves, and we wish to show in the paper that 
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this quantity is important for spontaneous emission. 
Therefore, from now on we assume fc 2 £ R. 
Modes /fe e are orthonormal, and their orthonormaliza- 
tion rules are: 

[d 3 re(r)f* x (r)-fUr)=Sex5(k-q), (8a) 

for radiation and substrate modes, while for guided 
modes: 

Jd 3 re(r)r qbX (r)-f kaC (r) = 

5 t \8{k x - q x ) 8{k y - q y ) 8 ab . (8b) 

In the considered case, relative electric permittivity 
e(r) = n 2 (r). Modes fke can be used to represent the 
generalized transverse delta function S^ 3 X (r,r'), which 
is a tensor conserving transverse fields. It is denned by 
the relation [7] 

J dVe(r') &(r,r') f ke (r') = fL(r) (9) 

(in this, and the remaining expressions in the paper, we 
use the summation convention for upper indices). The 
generalized transverse delta can then be represented 
as: 

%±(ry)=i fUr)fi:{r'). (10) 

Jk,c 

This function has two important properties, it is real: 



refractive index rii. Optical widths of these layers are 
defined as Ai = n\L\, A2 = niLi, in general, for the 
jth layer hS 3 ^ — ny)L z . Optical width of the whole 
elementary cell is then A = Ai + A2. It is known, that 
properties of photonic crystals depend on proportions 
of the structure, thus, to investigate its properties it is 
convenient to use parameters which are dimensionless, 
and we use A as a unit of length in the construction 
of such set. To describe the structure, we use refrac- 
tive indices, normalized optical widths A; /A (in gen- 
eral A"' /A), normalized position in the j = layer 
C = n (0) 2: /A' - 1 and normalized frequency = f/fo, 
where fa = c/A. To obtain more compact expressions, 
instead of the angle of incidence 1? we use its cosine 
r\ — cosi9. All expressions of our interest can be writ- 
ten in terms of these parameters, e.g.: 



k z z = 2-KrjOC,- 



A (o) 



±2jt0, 1 



V 2 A U) 



1 + 



U) 7 (o) 



- 1 



(14a) 
(14b) 

(14c) 



3. Hamiltonian of the system 

Electromagnetic field and a quantum system (e.g. an 
atom) interacting with each other can be described 
with the hamiltonian consisting of three parts: 



5l 3 *{r,r') = 5t 3 ,(r,r') 



(11) 



H = H CIn + Hint + H a 



(15) 



and 



(12) 



These identities can be easily proven using the relation 

/fee = -7e/-fee, with 



1. 



7e = 



for e = TE, 
-1, for e = TM, 



(13) 



which is satisfied by all modes, and Sfc 6 /fc tE = fk e , 
satisfied by substrate and guided modes, for which co- 
efficient Skc = rn exp(2ifc z L 2 ) [31]. 
Elementary cell of the considered photonic crystal con- 
sists of two layers, the first one with width L\ and re- 
fractive index m, and the second, with width L2 and 



Hem is the hamiltonian of electromagnetic field in in- 
homogeneous dielectric [7]: 



H- 



(16) 



where colons denote normal ordering, E and A are 
electric field and vector potential operators. H at is 
the hamiltonian of the quantum system. Because 
we are interested in the influence of the structure, 
we use the simplest possible expression, which for 
an electronic system, an atom or a molecule, in the 
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infinite nuclear mass approximation (so-called Born- 
Oppenheimer approximation for molecules) is the non- 
relativistic hamiltonian 

a 

where m is the electron's mass, p a is the ath electron's 
momentum operator and V is the potential operator, 
describing interactions of electrons and atomic nuclei. 
Formally, if denote eigenstates of H at with ener- 
gies Ei, it can be written 

F at =^S i |* i )(* i |. (18) 

i 

Hint is the interaction hamiltonian, which can be de- 
rived from the hamiltonian H&t through the minimal 
coupling p a ita = Pa — eA a , with e < being the 
electron's charge and A a = A(r a ) denoting vector po- 
tential in the position of the ath electron. In vacuum 

= Y,^ + V- ( 19 ) 

a 

The term emerging from is\ does not contribute 
to spontaneous emission, at least directly, because it 
contains no atomic operators and we neglect this term 
(what is in fact a usual procedure). Because = 
ihSabS 13 , momentum p a can be expressed as 

pj = g[ri,^]. (20) 

We wish to concentrate on interaction of electrons lo- 
calized in a small volume near an atomic nucleus (or 
nuclei in case of a small molecule), what is the case for 
e.g. dopant atoms in solid state lasers, where emission 
is related to transition between two atomic levels. It is 
possible to study transitions of non-localized electrons, 
like in conduction band of bulk semiconductor (see e.g. 
[34]), but that is not the case we are going to take into 
consideration in this paper. In the case of our inter- 
est, because optical wavelengths, which are of concern 
in case of light emitting devices, are much longer than 



the radius of an atom, the interaction Hamiltonian can 
be simplified by application of the dipole approxima- 
tion: A a ~ A(ro), where rn is the position of the 
system (understood as e.g. position of the center of 
mass). Then 

6£? = j:A t (r )[d i ,&«\, (21) 

where electric dipole moment d = ~}2 a er a . In dielec- 
tric, the dipole moment interacts not with the macro- 
scopic field, which is found as solution of Maxwell equa- 
tion, but with microscopic local field [7]. This effect 
can be accounted for by inclusion of field enhancement 
factor £dici, usually defined as the ratio of local and 
macroscopic electric field, but since the relation be- 
tween electric field and vector potential is linear, it 
can be used with vector potential as well. Thus, the 
interaction hamiltonian in dielectric medium becomes 

tfint = l -C iM A l {r ) [d ! , # at ] . (22) 

At this moment it is not necessary to define Ca\e\ ex- 
plicitly by an expression, and because its value depends 
on applied theory, we proceed further using the general 
symbol. The interaction hamiltonian obtained in the 
presented derivation is similar to the E ■ d form, but 
not equivalent. This form could be reconstructed if 
d, -ffmtj = 0, because the commutator in (22) would 
denote the time derivative, which could be transferred 
onto the vector potential, turning it into electric field. 
However, it is easy to check, that with the dipole ap- 
proximation \&, -ffmtj ~ — ifi-^-£diei-<4.(f o), where N 
stands for the number of electrons. Hence, we use the 
expression (22) for the interaction hamiltonian. 
In dielectric, vector potential and electric field opera- 
tors obey the commutation rule [7]: 

[A\r),Ei(r')] =_^ ± (r,r') (23) 

and can be in the usual way expanded into the effective 
resonator modes: 
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A(r)= ^ fc \[^- k ( akefUr) + a ' kJUr) )' (24a) 

E ^ =i ^ fe \[W ( ak ^W " , (24b) 

where afe e and a£ E are the photon annihilation and creation operators, obeying the bosonic commutation rules: 

[a k „a qX \ = [a^,aJ A ] =0, (25a) 
^afc t ,aJ A j = 5 c \5(k — q) (radiation and substrate modes), (25b) 
[afc a e, a^x] = 5t\8{k x — q x ) 8(k y — q y ) 5 ab (guided modes). (25c) 

Hamiltonian of electromagnetic field obtains then the form characteristic to an ensemble of harmonic oscillators: 

Hcm = Y^ huj k N kt , (26) 

J~k,e 

where Nke = a\ t a kt is the photon number operator with eigenvectors |. . .n ki . . .), where n ki = 0, 1,2, . . . and: 

ake | • • ■ n ke ■ ■ •) = y/nkl I ■ • ■ (rike - 1) . . .) , (27a) 
4 e \...n k€ ...) = Vn ke + l\...(n kt + 1) . . .) . (27b) 



4. Spontaneous emission 

Spontaneous emission is a process in which a relaxing 
atom emits a single photon. Thus, for the description 
of this phenomenon, the Hilbert space of ff a t can be 
restricted to only two eigenstates |*o) and (we 
assume that E\ > Eo), and the Hilbert space of if em 
to eigenstates |0) - with no photons, and |lfc e ) - with 
one photon in the mode with wave vector k and polar- 
ization e. For the whole system, the excited state 
is accompanied by the state with no photons |0) and 
the lower state l^o) - by field with one photon. The 
proper state vector is then: 

|V(*)) = C (t)|*i)|0)+^ C ke (t) |*o> \lke) ■ (28) 

The initial conditions adequate for spontaneous emis- 
sion are: 



In the two-level case, states |\l/o) and can be 

chosen so that the transition dipole moment dio = 
(*i|d|*o) £ R 3 . In this paper we calculate the 
spontaneous emission rate following the procedure de- 
scribed in [35]. After substitutions: 

C (t) = b (t)e X p^- l -^j , (30a) 
C fee (t) = b kt (t) exp (-i ^ + ^ t) , (30b) 

where bo(t) and b kt (i) are the probability amplitudes 
in the interaction picture, the equations of motion be- 
come: 

9h M)e i(ni, ^ 1 (31a) 
^| i =-.9L6o(i)e i( "'=- nio)t , (31b) 



Co(t) = l,C ke (t) = 0. (29) 



where angular frequency of the atomic transition fiio = 
(Ei - Eo) /h and 
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In this case, if ro denotes the position of the quantum Formal solution of (31b) is: 
system: 

^dicifiiudin • fke(ro) ,„„.. r* ., _ . 

9kt = • (M) bk '® = ~ 9L / dr6o(r)e^- n "^, (34) 

j o 

Symmetry of multilayer structure allows to put ro = 

^ g_ therefore: 

I 

r ) e i(n 10 -u, fc )(t-r) > (35) 



Assuming that 6n varies slowly during the time interval from to t and pulling it in front of the time integral 
(Markov process approximation) one obtains: 

1T = -^ \9k.\'bo(t)J*dr^"- u ^\ (36) 



In the limit t — > oo the integral over r becomes: 



lim / 

t->ooy 



^giCnio-WfcJCt— r) _ 



7T<5(fiin — Wfe) - P 



filfj — OJfc 



(37) 



where V denotes principal value. The term containing 
the delta function describes the decay process and can 
be used to determine the spontaneous emission rate, 
while the principal value term results in line shift 5uj, 
which is not of our concern here. In the long time limit 



+ i<5w b (t), 



therefore 



where 



dbo 
IF 



bo(t) - exp( -y ) e lSu 



F = Trm + Tsm + Tga 



(38) 



(39) 



(40) 



is the spontaneous emission rate. Contributions from 
different kinds of modes to F are: 

r R ,M=27rV/ d 3 fc \g kc \ 2 S{n 10 - to k ) , (41a) 

£ J RM 

r SM =2TvJ2 [ d3fc \9k,\ 2 s{n w -tu k ), (4ib) 

. JSM 



F GM =2tt^ ^d 2 fc|| ^ 

e aGGM(*!||,e) 
X \dk a e\ 2 6(0.10 — Wfc a ) ■ 



(41c) 



In case of free space, expression for F leads to the well- 
known Weisskopf-Wigner rate: 



3nhc 



(42) 



With the knowledge of fen, solution of the equation of 
motion for b kt is: 



„ 1 - exp(i (Suj + o; fc - n 10 ) t - ^) 
bke(t)=g kt — ^ — ■ — - p ^ , (43) 



i (Su) + uj k — Oio) — § 



thus, in the limit t — ► oo: 



(44) 



5. Numerical calculations 

The integration of \g k c\ 2 necessary to calculate numeri- 
cal value of F is best performed in spherical coordinates 
(k, $, ifi), where is the angle of incidence and ip is the 
azimuthal angle. Results of integration over p can be 
obtained easily. Every /cte lies within the xy plane 
and can be written in function of p: 

/fcTE = -Pke x simp + P k e y cosip, (45) 



where 



Pk — e y ■ /fcTE 



(46) 
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Because the structure has rotational symmetry, axes x Then: 
and y can be chosen in the way that 



Then: 



dio = d«e x + d z e z 



r2w 

/ dp |dio • /fcTE] 2 = vrdj| |P fe | 2 . 
Jo 



(47) 



dip | dio • fk 



ird x \Q k »\ +27rd z | 



(51) 



Integration over $ has to be performed numerically. 

5.1. Radiation modes 

(48) 

Radiation modes spread in the range of < i? < #£ 
and 7r — &s < $ < i", where 



Every /i-tm lies in the plane of incidence and can be 
written in function of ip: 

Atm = Q k \\e x cosp + Q kl] e y smp + Q k± e z , (49) 

with 



fs = arcsm ■ 



{ n (Q), n (L), n(R)} 
n (0) 



(52) 



Contribution of radiation modes Trm to spontaneous 
emission rate is: 



Qk\\ — e x ■ fk 



I <p=0 



and 



e z • fkTM- 



(50a) 

(50b) 
I 



27rn( 0) f2 2 



V / dtf sim? / d^ [\g ke 
~ Jo Jo 



fc = Tl( )Q 10 /c ' 



(53) 



In the second integral over ip it is possible to change the variable to ip + n and obtain |<7-fce| instead of |gfe» e |. 
Then, because |<7-ft e | = \gke\, using the known results of integration over ip, after a few operations one arrives at 



d? 



d v [d'i(\P k \ 2 + \Q kll \ 2 )+2d 2 z \Q k± \ 2 ] ( 



with normalized transition frequency: 



27T/o 



rjs = cosi?s and 



|Pfe|2 = PTE(fe) (i+Be^^^J) 

n (0) V 1 J/ 



\Q k± f = (!_„») El^l (l + Re{cTM(fc)e 2ifc -«}) 



For radiation modes 



(l-\r L \ 2 )e- 2ik * L * +r L jl-\r R \ 2 ) 
l-\r L r R \ 2 



(54) 
(55) 

(56a) 
(56b) 
(56c) 

(57) 



Coefficients (56) can be easily calculated using explicit definitions of field distributions from [29] or [31] (we do 
not recall them here, because it would be necessary to refer to a significant part of calculations from [29] in order 
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to explain the symbols used in them), knowing, that in 
significantly simplified, i.e.: 



those papers a few formulas for radiation modes can be 
I 



F kf = 



"(0) 



'(()) 



1- \R 



\Uk)Y 



(58a) 
(58b) 

(58c) 



In free space, where all the modes are of radiation type, 
(54) correctly reproduces the Weisskopf-Wigner spon- 
taneous emission rate, and in a multilayer structure, 
contribution of each radiation mode is proportional to 
its mode spectrum. 

5.2. Substrate modes 

Substrate modes are found in the ranges of #s < ■& < 
•&G and 7r — i?g < •& < n — i?s, where 



^ = arcsin min ^ (0) ' maX ^ W ' n(fl) ^. (59) 

n (0) 



However, these ranges contain the same modes [31], 
therefore they must not both enter into the integration. 
Thus, the contribution to spontaneous emission rate is 
given by: 



Tsm = ■ 



d 10 



I" 

■J lir 



dr/ 



X 



[dl(\P k \ 2 + 



+ 2d 2 z \Q k± \ 2 



e=e 10 



where t)g = cosQg- For substrate modes: 
id ,2 „1 + Re{f T E(fc)exp(2ifc 2 2 )} 

\^k\ = 2 = ' 

, ,2 _ 2 1 - Re{gTM(fc) cxp(2ifc z z )} 

-TfcTM 

IQfcxI 2 = 2 (1 - T7 2 ) 

1 + Re{gTM(fc) exp(2ifc z z )} 

-FfcTM 



(61a) 
(61b) 

(61c) 



where £ e (fc) = r* R exp(— 2ik z L z ), with the assumption 
that < ?i(L), and the coefficient F kt is given by 

(58a), just like in case of radiation modes, thus, these 
expressions are also proportional to mode spectrum. 



5.3. Guided modes 

Guided modes occur at discrete angles # from the 
range #g < "& a < tt/2 or 7r/2 < $ n < n — $g. Just 
like in the case of substrate modes, guided modes with 
•da and n — -d a (and the same polarization) are in fact 
the same mode, therefore summation has to restrict to 
only one of these angles. Angles of incidence of guided 
modes are determined by the relation 



i<h a 2ifc z 



(62) 



which allows to find them in function of frequency, or 
alternatively, their cosines in function of normalized 
frequency rj a (0). Because 



2im m 0ri a (0) 
A ' 



(63) 



for a guided mode |7\l| = |rn| = 1 and 



d& 



(64) 



it is easy to find, that 



(new dr,a 

(60) de=~ 



Im{ 7I ^^(r L r«)}+4^41 



Illl 



{^ v (r L r R )} 



+ 4tt0 



A(°> 



(65) 



This equation can be solved numerically and allows 
to obtain values of r) a and its derivative for different 
easier, than by solving (62) directly (with e.g. bisection 
method). 

Because 

2nn m 0^1-tf(0) 



A 



(66) 



integration over fcy can be replaced by integration over 
0, which is carried out immediately, as 



5(Qio — uj ka ) = 



5(0 - gi 
2vr/ 



(67) 



and the guided modes' contribution to spontaneous 
emission rate turns out to be: 
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37r 2 n( )£d iel r fs 
2&To 



E 



aSGM QiTE 



4a|p,| 2 

df,, 



1 - I) o ,te(0) ~ Orj a ,TE{0) 



+ E 



d 2 A|Q fc |||- + 2d 2 A|Q fcJ 
d? 



d<9 



i - 1b,tm(0) - 9l)«,TM(e) 



rf??a,TM(@) 
dO 



© = ©10 



where GMe.e denotes the set of guided modes for nor- 
malized frequency O and polarization e, For guided 
modes Q^y and Q^x are given by expressions of 
the same form as in the case of substrate modes and 
the relevant definition of can be found in [30] or 
[31] as well. However, it is worth to stress, that in 
case of guided modes coefficients Fke have the dimen- 
sion of length, thus ratios Fke /A are dimensionless. 



dVaA&X\ 



is often ne- 



The term 1 - j?L(<9) - Or) a , e {&) 
glected, e.g. in [21]. It is in fact equal to one in a res- 
onator with metallic mirrors, but in a multilayer struc- 
ture it can differ from unity quite significantly. 

6. Exemplary results 

The expression derived in previous sections allow to 
calculate and present numerical values of the sponta- 
neous emission rate, provided an explicit expression for 
£dici is specified. Following Glauber and Lewenstein 
[7], we put 



■Cdiel — 



3n 



(») 



2n 2 0) + l' 



(69) 



what is the proper choice for substitutional impurities 
[36]. We have conducted calculations for the following 
structures: 
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Figure 2. Absolute value of reflection coefficients, mode spec- 
trum (both for rj — 1) and spontaneous emission 
rates for structure SI. 



51, with Ai/A = 0.65, ni = 3, ni = 1.6, n^ R ) — 
ri(L) = 1 and Nr = — Nl = 12 (i.e. 6 periods in 
each direction starting the count from the j = 
layer); 

52, with Ai/A = 0.5, m = 1.4, n 2 = 1.2, n {a) = 
1, n (L) = 1.1 and Nr = -N L = 30; 

53, with Ai/A = 0.5, width defect in the j = 
layer: A (0) /A = 0.7, m = 1.4, n 2 = 1.2, n (R) = 
1, n (L) = 1.1 and Nr = -Nl = 20 (i.e. 10 
periods in each direction starting the count from 
the j = layer). 



Results for structure SI are presented in Fig. 2. The 
top plot contains the absolute values of reflection co- 
efficients of stacks of layers on both sides of the j = 
layer, calculated for rj = 1, i.e. perpendicular incidence. 
The middle plot contains mode spectrum (normalized 
to the free space value) for the same angle of incidence. 
These two plots show the band gaps forming in the 
photonic crystal. Band gaps appear near the normal- 
ized frequency 0b = 0.5, corresponding to the Bragg 
wavelength Ab = 2A and its multiplicities. For a band 
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Figure 3. Absolute value of reflection coefficients, mode spec- 
trum (both for 77 — 1) and spontaneous emission 
rates for structure S2. 



Figure 4. Absolute value of reflection coefficients, mode spec- 
trum (both for r\ — 1) and spontaneous emission 
rates for structure S3. 



gap to form, a high quality factor of the cavity (layer) 
is required, therefore they appear in the regions with 
high reflection, which are characterized by low values 
of mode spectrum. In perfect photonic crystal, mode 
spectrum in a band gap would be equal to zero, but in 
a real, finite structure it is not possible. Mode spec- 
trum is high at the edges of a band gap, resembling 
an effect of the modes being "pushed out" of the bang 
gap. It is also possible to relate values of mode spec- 
trum to profiles of field distribution - for high value, 
the mode has strong field in the layer, while for low val- 
ues its field in the layer is weak. Because the emission 
rate is high for modes exhibiting strong fields at the 
location of the atom, modes with high mode spectrum 
are expected to share the most of emitted power, while 
those with low mode spectrum have negligible contri- 
butions to the overall rate. These observations allow to 
interpret the bottom plot, showing emission rate split 
into contributions from all types of modes. Because for 
structure SI the dielectric outside is the same, there are 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Figure 5. Contribution to T at 77 = 1 for structure S2. 



no substrate modes and Fsm = 0. The band gap, seen 
in the middle plot, for a lower rj would slide towards 
higher (because the k z component of the wave vector 
is determined by Ab), but in this structure band gaps 
for all r) partially overlap, thus, it exhibits a full band 
gap. That is why in the plot of the radiation modes' 
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Figure 6. Contribution to T at J7 = 1 for structure S3. 



contribution Frm/F there can be seen a region of <9, 
in which the emission is practically forbidden. How- 
ever, there is present a contribution from a few guided 
modes, meaning that for transition frequencies in the 
band gap the excitation of an atom would decay quite 
slowly with emission into these modes. 
Similar results for structures S2 and S3 are plotted in 
Figs. 3 and 4. For these structures, m and ri2 are 
too low for a full band gap to form, therefore they 
do not inhibit spontaneous emission as selectively as 
SI. Instead of trying to trace characteristic spots in 
r, which is roughly constant, it is much more interest- 
ing to investigate the contributions to the decay rate 
from modes with a given value of 77 = 1 (the same 
as for which mode spectrum was plotted), shown in 
Figs. 5 and 6. These plots make it evident, that there 
is a strict relation between mode spectrum and contri- 
bution to decay rate, and therefore probability of emis- 
sion into a particular mode, given by (44), as well. It is 
also interesting to note, that the defect of structure S3, 
which results in a defect mode seen in the mode spec- 
trum inside the band gap, does not increase the total 
decay rate, but it rather helps to direct the emission 
into specific modes. Because at each frequency defect 
modes occur only at one angle of incidence, while max- 
ima at the edges of band gap can appear at two, thus, 
structure S3 makes a better control of the emission, 
than S2. 

Characteristic emission rate for an atom in a specific 
layer is proportional to mode spectrum in that layer, 
therefore analysis of variation of mode spectrum in dif- 
ferent layers reveals important information about the 




Figure 7. Mode spectrum of a defect mode in consecutive n\ 
layers of photonic crystal SI with a defect intro- 
duced in layer jV dc f (for rj = 1, N stands for the 
index of elementary cell). 



structure. One of interesting results is presented in 
Fig. 7. In this figure, there is plotted mode spectrum 
of a defect mode in n\ layers of different elementary 
cells of structure SI with a defect introduced in one of 
its layers. It is clearly seen, that the defected layer is 
characterized by the maximal value of mode spectrum, 
which becomes quickly much smaller in surrounding 
layers. The maximal value depends also on the situ- 
ation of the defected layer. It follows from obtained 
expressions, that the same is the behaviour of sponta- 
neous emission rate into the defect mode. Similarly, 
in case of a ID photonic crystal with multiple defects, 
analyzed in [32], the results presented therein can be 
immediately related with spontaneous emission rate. 
This conclusion could not have been drawn from den- 
sity of states calculated for the whole structure or with 
a model assuming external excitation and is a good ex- 
ample of one of the effective resonator model's advan- 
tages. 

7. Summary 

In this paper we have obtained the expressions de- 
scribing spontaneous emission rates in a structure of 
finite one-dimensional photonic crystal with arbitrary 
defects. Our derivation has been based on the effective 
resonator model, which allows to calculate field distri- 
butions and defines a quantity called mode spectrum, 
which contains information about physical properties 
of a particular layer of the structure. Thanks to this 
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property mode spectrum is much more useful for anal- 
ysis of the structure's properties, than e.g. density of 
states, characterizing the whole structure, because it 
allows to easily observe effects associated with differ- 
ent parts of the structure. We have defined a set of 
dimensionless parameters of the structure, which are 
very convenient to work with, in particular the derived 
formulas for the spontaneous emission decay rate are 
fully expressed by them. We have discussed contribu- 
tions from modes of different types and shown, that 
contributions from radiation and substrate modes are 
proportional to mode spectrum. This is an important 
result, because calculation of mode spectrum is much 
easier than of field distribution and provides an easy 
way of investigation how a structure affects emission of 
electromagnetic radiation. Employing mode spectrum 
for an analysis is particularly beneficial in case of de- 
fected structures, allowing to easily reveal defect modes 
or their behaviour throughout the structure, what is 
not always easy to accomplish with other models. Ex- 
emplary results included in the text indicate, that the 
developed description can easily characterize sponta- 
neous emission from defected layer or any nearby layer, 
and provides information on how localization of field 
changes when the defected layer is moved in the pho- 
tonic crystal. Thus, we have shown that the effective 
resonator model is a tool suitable for modeling of spon- 
taneous emission from multilayer structures, includ- 
ing defected finite one-dimensional photonic crystals, 
which can be easily applied in various practical designs. 
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